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ABSTRACT 



A new generation of low frequency radio telescopes is seeking to observe the 
redshifted 21 cm signal from the Epoch of Reionization (EoR), requiring innova- 
tive methods of calibration and imaging to overcome the difficulties of widefield 
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low frequency radio interferometry. Precise calibration will be required to sepa- 
rate the small expected EoR signal from the strong foreground emission at the 
frequencies of interest between 80 and 300 MHz. The Moon may be useful as a 
calibration source for detection of the EoR signature, as it should have a smooth 
and predictable thermal spectrum across the frequency band of interest. Initial 
observations of the Moon with the Murchison Widefield Array 32 tile prototype 
show that the Moon does exhibit a similar trend to that expected for a cool 
thermally emitting body in the observed frequency range, but that the spectrum 
is corrupted by reflected radio emission from Earth. In particular, there is an 
abrupt increase in the observed flux density of the Moon within the internation- 
ally recognised Frequency Modulated (FM) radio band. The observations have 
implications for future low frequency surveys and EoR detection experiments 
that will need to take this reflected emission from the Moon into account. The 
results also allow us to estimate the equivalent isotropic power emitted by the 
Earth in the FM band and to determine how bright the Earth might appear at 
metre wavelengths to an observer beyond our own solar system. 

Subject headings: dark ages, reionization, first stars — Extraterrestrial intelli- 
gence — Moon — Techniques: interferometric 



1. Introduction 



The astrophysics of the Intergalactic Medium over the redshift interval z f» 200 to 8 is 
expected to give rise to observable p h enomena in the reds h ifted 21 cm line of neutra l hydr ogen 



^see e.g. IPritchard & Loebl (1201 if ): Morales fc Wvithd (120101 1: iFurlanetto et all (120061 ) for 
detailed reviews). The 21 cm line from this period, which extends from the Dark Ages 
through the Epoch of Reionization (EoR), redshifts to radio frequencies of between 30 and 
200 MHz, making it a suitable target for the next generation of low frequency radio telescopes. 
Use of the redshifted 21cm line of neutral hydrogen to study the EoR will require exquisite 
calibration in order to achieve milli-Kelvin level spectral line sensitivity in the presence of 
Galactic and extragalactic foregrounds that are brighter than the expected signal by at leas t 
2 to 3 orders of magnitude dshaver et al 1999 : Rogers fc Bowman 2008 : Bernardi et al 2009 ). 



One suggestion (Sh aver et al 



1999) for in strumental calibration for observation of the global 



EoR signature ( jBowman fc Rogers! 120101 ) would be to use the smooth, thermal spectrum of 
the Moon for comparison with the sky brightness. 



The goals of the work reported here are to assess the suitability of the Moon as a 
comparison source for EoR detection and to test novel observing modes for the Murchison 
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Widefield Array (MWA) (ITingay et alll2012. in presa ; lLonsdale et alll2009f ) that would span 
the full range of the uncertainty in redshift that surrounds the reionization phase change. 

A concern at the outset has been that the lunar spectrum could be corrupted by reflected 
communications and navigation trans missions from the Earth. It has been known since the 
US Army Signal Corps' Project Diana (IDeWitt fc Stodolalll949l ) that radio signals could be 
reflected by the Moon and received on Earth. Indeed, it is now com mon practice for amat eur 
radio enthusiasts to bounce communications signals off the moon ( iKatz fc Franco! l201l[ ) to 
be received by listeners on t he other side of the planet. Radio leakage from the Earth has 
been modelled theoretically (I Sullivan et all 119781 ) and measured by observin g the reflected 
signals off the Moon using the Arecibo telescope (jSullivan &: Knowleslll985f ). However, the 
use of radio communications on Earth has changed since these studies were undertaken and 
the EoR waveband of 80 to 200 MHz covered by the MWA deserves an updated study 
and continued monitoring. A strength of the present observations with the MWA 32 tile 
prototype (MWA 32T) is that the gain of the synthesised beam isolates the lunar reflections 
from any ambient signals that may enter weakly into the radio quiet zone of the Murchison 
Radio Observatory (MRO) on more direct paths from local transmitters. 

At the frequencies of interest the radio sky is bright and dominated by synchrotron 
radiation f rom the Milky Way. The sky background varies as a power law function of 
frequency ([Rogers fe Bowman! 12008c Ide Oliveira-Costa et all 120081 ) and the sky brightness 
temperature is given by: 



T, 



sky 



L . r )0 



V 



150 MHz 



K 



wher e the spectral index at high galactic latitudes has been measured ( IRogers fc Bowman 
20081 ) to be 2.5 and at the sky position of th ese observations (J2000 coordin ates h 51 m 6 s , 
11°24'40") the sky temperature T 150 is 248 K jLandecker fc Wielebinskilll97ol ). 



The brightness temperature of the Moon can be modelled by a comb ination of three 
facto rs; thermal emission, assuming a Moon brightness temperature of 230 K (jHeiles fc Drake 
1963), reflected radio sky background emission, assuming a Moon reflectivity of 7% ([Evans 
19691 ) and sky temperature given by equation (pQ), and reflected emission from the Earth. 
Reflected e mission from the Sun , even at Full Moon, is expected to contribute less than IK for 
a quiet Sun (IHagfors et allll969l ). This reflected solar emission is not included since, based on 
the information available from the National Oceanic and Atmospheric Administration Space 
Weather Prediction Center (NOAA SWPC), the solar activity during our observing period 
was low. The contribution to the brightness temperature from reflected emission from the 
Earth is left as an unknown parameter in this model as it is not well understood at these 
frequencies and for recent times. These observations aim to contribute to the characterisation 
of this reflected emission. 
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Brightness temperatures (T& in K) are converted to flux densities (S in Jy) by: 

where k is the Boltzmann constant in J K _1 , Q is the solid angle subtended on the sky in 
steradians and A is wavelength in m. The value of Q used is 6.0 x 10 -5 which corresponds 
to the Moon's angular diameter of 30.0 arcminutes at the time of these observations. The 
expected difference in flux density between the Moon and the sky, ignoring reflected emission 
from Earth, is plotted in Figure [1] and shows that at the lower end of the frequency band the 
Moon should appear as a 'hole in the sky', crossing over to being a positive source relative 
to the sky background at around 150 MHz. Recovering this Moon-sky difference signal was 
an important aim of the MWA observations described in this paper. 



l — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — r 




I i i i I i i i I i i i I i i i l_ 

80 100 120 140 160 

Frequency (MHz) 



Fig. 1. — Theoretical difference between Moon and sky flux density at J2000 coordinates 
h 51 m 6 s , 11°24'40", for observations at 0030-0132 LST on UTC dates 24 and 25 September 
2010 respectively. 
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Observations 



The MWA is a new type of radio interferometer array that is designed to observe the 
redshifted 21 cm sig nal at frequencies of between 80 and 300 MHz, corr esponding to a redshift 



range of z=4 to 17 ( ITingay et all 12012. in press! ; lLonsdale et alll2009l ). It consists of a large 



number of small, crossed-dipole antennas grouped into 'tiles' whose beams are electronically 
steered using analog beamformers. Each tile is made up of 16 dipole antennas arranged 
in a 4 by 4 grid above a conducting mesh ground plane. A 32 tile prototype instrument 
(MWA 32T) was operated over the period September 2007 to September 2011 within the 



and lWilliams et all (120121 ) 



Oberoi et al 


( 


2011 


), 


Ord et al 


( 


2010) 



The observations were conducted with the MWA 32T during Expedition 14 to the MWA 
site, which took place in September 2010. The 32 tiles of the prototype system were arranged 
as shown in Figure El The positions are relative to the array centre. This tile configuration 
gives an instantaneous UV coverage as shown in Figure [3] for a pointing at the Moon at 
azimuth 357 degrees, elevation 52 degrees and frequency 80 MHz. Multi-frequency and 
Earth-rotation synthesis were used to maximise UV coverage for imaging. With this tile 
configuration the maximum baseline is approximately 350 m, giving a spatial resolution of 
approximately 30 arcminutes at 100 MHz. 

The Moon observations were made in a novel 'picket fence mode' which reallocated the 
30.72 MHz of signal processing capacity to sample the frequency range from 80 to 200 MHz 
with 24 evenly spaced 1.28 MHz wide sub-bands. Four different sub-band configurations 
were used to fill the whole spectrum across this frequency range. The 24 sub-bands of the 
MWA 32T are further divided into 32 fine channels, giving a spectral resolution of 40 kHz. 

Observations of the Moon were made over the identical range of 0030-0132 LST on UTC 
dates 24 and 25 September 2010 (approximately 1630-1732 UTC). The resulting dataset con- 
sisted of 32, 107 second integrations (eight scans for each channel configuration) on both days. 
The telescope was configured such that the beam centre was aimed at the computed position 
of the Moon at 0100 LST on 24 September (J2000 coordinates h 51 m 6 s , ll o 24'40") and the 
field allowed to drift across this fixed beam position on both days. Identical beamformer 
delay settings were used on both days. This was done to maximise the stability of the system 
as re-pointing the beam has complicated gain effects that are yet to be fully characterised. 

A calibration dataset, using a setup identical to the Moon scans, but with the fixed 
telescope pointing centred on Fornax A, was also taken in the hour following the Moon scans 
on each day. These data were used for obtaining initial phase calibration solutions in order to 
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Fig. 2. — Positions of MWA 32T tiles relative to the array centre. 
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Fig. 3. — MWA 32T instantaneous UV coverage for a pointing at the Moon at azimuth 357 
degrees, elevation 52 degrees and frequency 80 MHz (wavelength 3.75 metres). 
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start the self-calibration pr ocess. Fornax A could not be used to set the overall flux density 
scale as it has been found f jOrd et a]||2010l ) that while phase solutions transfer well between 
tile pointings, the amplitudes of the complex tile gains cannot be applied to different tile 
pointings without being in error by several percent. Fortunately, the relatively bright source 
3C33 lies in the main beam close to the Moon's position and this source was used to tie the 
Moon measurements to an established flux density scale. 



3. Calibration and Imaging 

The raw correlator data were averaged in time by 6 seconds and converted to UVFITS 
format. The dataset was then imported to the Astronomical Image Processing System 
(AIPS) for data reduction. The UV data were examined at the finest spectral resolution 
and bad channels flagged manually. Radio Frequency Interference (RFI) was sparse and 
generally limited to a few discrete narrow bands. Four fine channels at the edge of each 
sub-band, as well as the central fine channel were also flagged as these were known to be 
corrupted by aliasing and rounding error, respectively. Tiles 2 and 6 were also flagged during 
the observations as they were not behaving as expected. In order to reduce the subsequent 
processing time, the UV data were then averaged in frequency to give 0.64 MHz sub-bands. 
The complete dataset, including both days, was concatenated into one AIPS file to allow 
self-calibration of the dataset using the same sky model for both days. This approach tied 
the tile gain solutions for both days to a common reference, which was important for the 
Moon-sky differencing experiment. 

A subset of the Fornax A data were used as an input to the AIPS 'CALIB' task to get 
an initial set of complex tile gain solutions for calibration. A calibration was performed on 
these data, assuming a sky model of a point source located at the centre. The UV range 
in the initial calibration was selected to exclude longer baselines so that Fornax A would 
not be resolved. An image was then made using the AIPS task TMAGR', im plementing 



Clark's variation of the Hogbom CLEAN algorithm ( IClarkj|l980[ ISchwabl Il984l ). Faceting 
was used to account for wide-field effects. The clean component list generated by this imaging 
process was then fed back into the calibration procedure as an improved sky model. This 
self- calibration process was repeated for several iterations using the full UV range and phase 
solutions only, until convergence was achieved. A final phase and amplitude self-calibration 
iteration was performed and the resulting image of Fornax A clearly showed a bright double 
source, with many fainter radio sources in the field of view. 

The tile gain solutions from the final self-calibration iteration of Fornax A were then 
applied to the Moon data as a starting point for self-calibration. The same procedure as 
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described above for Fornax A was used to generate an image cube using the combined set of 
Moon data from 24 and 25 September. Imaging was performed in 10 MHz bands to maintain 
a high signal to noise ratio, resulting in an image cube with 9 frequency slices, centred on 
the Moon's position at 0100 LST on 24 September. After the final calibration iteration 
was complete, the tile gain solutions were applied and the data for both days were imaged 
separately, producing an image cube for each day. Using the same sky model to self-calibrate 
the data for both days at the beginning of the process effectively tied the tile gain solutions 
to a common reference, allowing a fair comparison of the two days when imaged separately 
at the end. The motion of the Moon during the hour of observations on each day is 19 
arcminutes, which is approximately the size of the synthesised beam at 150 MHz (taking 
the diameter of the synthesised beam as A/350m). So at the higher end of the observed 
spectrum we would expect some smearing of the Moon image due to this motion. 

Examination of the cubes revealed that the Moon is present as a bright source in two of 
the image slices, corresponding to 94 MHz and 104 MHz, on both days. On 24 September 
(Figure HI left panel), the bright source appears at the centre of the field and on the 25 
September (Figure HI right panel), the bright source has moved approximately 13 degrees 
across the sky as expected of the Moon. The result shown in Figure H] where the Moon 
appears to be a bright radio source rather than a cool 'hole in the sky' does not agree with 
the simple prediction from Figure [U which neglects reflected emission from the Earth. 

To investigate this spectral behaviour in more detail, a higher spectral resolution image 
cube was generated for the entire 50 degree diameter field using 1.28 MHz channels. We 
refer to this data cube as the Intermediate Calibration Cube, since the relative flux densities 
for the sources at each frequency are well determined but the overall gain passband and flux 
density scale have not yet been imposed. The spectrum of the integrated flux density of 
the central 'Moon' pixels in the Intermediate Calibration Cube was plotted and is shown in 
Figure It is clear that the Moon's brightest disk averaged emis sion correspo nds exactly 



with the internationally recognised FM radio transmission band (jACMAl 120091 ) of 87.5 to 
108 MHz. 

The imaging described up until this point has used self-calibration with an arbitrary 
flux density scale. Setting the flux density scale for the MWA 32T is difficult due to the 
currently uncharacterised behaviour of the overall gain of the system for different beamformer 
settings. This means that the traditional method of pointing at a known flux calibrator and 
then pointing at the target source cannot be used reliably and a known source within the field 
of view of the target is required. The relatively bright source 3C33, located approximately 
4.7 degrees from the pointing centre, was therefore used as a calibration reference to establish 
a flux density scale for the Moon. 3C33 exhibits a power law spectrum with measured flux 
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Fig. 4. — Moon field images on 10 MHz bandwidth, centred at 94 MHz, for 24 September 
2010 (left panel) and 25 September 2010 (right panel). The flux scale is arbitrary at this 
stage of the calibration process and no primary beam correction has been applied. Darker 
greyscale indicates higher flux density. 
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Fig. 5. — Moon spectrum at 1.28 MHz resolution on 24 September 2010 with an arbitrary 
flux density scale. The vertical lines denote the FM radio band of 87.5 to 108 MHz. 



densities of 105.3 ± 0.4 Jy at 74 MHz flKassim et all 120071 ) and 58.3 ± 2.9 Jy at 178 MHz 
( jKiihr et all Il98lh . Interpolating these data using a power law with a single spectral index 
gives the following equation for the flux density of 3C33: 



S C alc{v) — S74 



74 MHz 



v 




(3) 



where v is the frequency in MHz, S74 = 105.3 ± 0.4 and a = 0.67 ± 0.06 

The spectrum of 3C33 extracted from the Intermediate Calibration Cube was measured 
using the AIPS task 'ISPEC. A correction was made to this measurement to account for 
the drop in gain of the primary beam as a function of frequency for the source which is 
approximately 4.7 degrees from the beam centre. This rudimentary primary beam correction 
was based on a simple analytical model that didn't take into account the different beam shape 
of the two dipole polarisations, but is accurate for positions close the pointing centre. We 
refer to this corrected flux density of 3C33 as Sicc{ v )- 

Figure [5] does not show the smooth spectral shape expected for a thermal spectrum. 
Even outside the FM band there is considerable structure above the nominal noise level. 
The flux density calibration specified by equation (J3j) is smooth on these scales and cannot 
cause this sort of structure. 

Since the observations were made on two days with identical tracks as a function of 
sidereal time, certain classes of systematic imaging artefacts are expected to cancel when 
the difference between the two days data is examined. Specifically, the chromatic effects 
of residuals left from bright sources and integrated source confusion, observed with the 
incomplete UV coverage attained with the MWA 32T, can couple to spectral structure that 
will vary from beam to beam throughout the image cube. As a test of the presence of 
artefacts, an identical imaging and calibration process was used to produce an ICC for 
the 25 September, and then a difference image cube was constructed by subtracting the 
25 September cube from the 24 September cube. Figure [6] shows an image slice from the 
difference cube centred at 99.835 MHz. 

For a perfect instrument and perfect calibration on both days of observation with identi- 
cal conditions you would expect the RMS noise in the difference image (Figure [6]) to approach 
the thermal noise level, since all sources (including low level confused sources within the syn- 
thesised beam) would be perfectly subtracted. 



The expected noise level in the Moon images, a in Jy, is given by: 



- 13 - 




02 00 



01 30 



00 00 30 

RIGHT ASCENSION (J2000) 



Fig. 6. — Image slice from the Moon-sky difference cube centred at 99.835 MHz and covering 
0.92 MHz of bandwidth (since 28% of the sub-band is flagged). The flux scale of the image is 
arbitrary at this stage of the calibration process. Darker greyscale indicates higher (positive) 
flux density. 
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(adapted from iThompson et all (120011 ) equation 6.62) where /cb is the Boltzmann constant 
in J/K, T = T&y + Tr rv 700 K at 100 MHz, A e g is the effective collecting area of a tile w 
21.5 m 2 ( Tingay et alll2012. in press! ). N is the number of tiles (N= 30 for these observations), 



e c is the correlator efficiency (assumed unity), B is the instantaneous bandwidth in Hz, n p 
is the number of polarisations and r is the integration time in s. Equation (jlj) assumes 
natural weighting of the visibilities and is a close approximation to the noise properties of 
the weighting scheme used by default in AIPS. 

We find the RMS in the difference image cube to be higher than the expected thermal 
noise given by equation (jlj) by a factor of between 4 and 7 across the frequency range 
and that the excess is highest in the FM band. For an imaging sub-band of 0.92 MHz 
(since 28% of fine channels are flagged) and an integration time of 7 minutes (the hour of 
observations was split into four channel configurations and half of the integration time was 
lost due to network bandwidth limitations) the RMS in the difference image slice shown in 
Figure [6] is approximately 1.3 Jy, whereas the expected thermal noise given by equation (jlj) 
is approximately 0.2 Jy (including an additional factor of \/2 due to the differencing). The 
main factor causing this discrepancy is thought to be imperfect calibration on both days 
causing incomplete subtraction of sources in the difference image. The same sky model was 
used to calibrate the concatenated data from both days. However within the FM band, 
there is a 60 Jy source (the Moon) in two different positions 13 degrees apart on subsequent 
days. This leads to a calibration solution that is slightly in error for both days and a higher 
than expected RMS noise in the difference cube. Another factor could be differences in the 
ionosphere between the two days which would affect the positions of sources in the field, 
resulting in imperfect subtraction. The errors in the difference image are included as error 
bars in Figures [8] and [TD1 

The integrated flux density of the central pixels of the difference cube gives the difference 
between the Moon flux density on the 24 September and the sky flux density on 25 September. 
This Moon-sky difference spectrum was calibrated by application of the correction factor 
Scaiciy)/ Sicciy) to give the calibrated spectrum shown in Figure These observations 
clearly show the Moon to be a source of approximately 60 Jy across the FM radio band from 
87.5 to 108 MHz and many of the residual spectral features apparent in Figure have indeed 
been removed by the differencing process. The uncertainty in the flux density of 3C33 leads 
to an uncertainty of less than 3% in the flux density calculated for the Moon. Following 
the same steps, we extracted a lower signal-to-noise ratio spectrum for the second day, when 
the Moon had moved approximately 13 degrees from the centre of the tile primary beam. 
Since the tile primary beam Full Width Half Maximum (FWHM) varies as approximately 
25°/(* / /i50MHz), the Moon's position on the second day corresponds approximately to the 
half power point of the tile primary beam at 150 MHz. This spectrum is consistent with the 
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spectrum obtained while the Moon was at the centre of the beam as shown in Figure [7J 

Excluding the FM radio band, the measured Moon-sky difference has the negative to 
positive trend of Figured], however the slope of the measured difference plot is still not a good 
fit to the predictions which have neglected reflected emission from Earth as shown in Figure 
[HJ The error bars in Figure [8] represent the RMS noise in each image slice for a square region 
of the difference image in between the Moon's position on each of the subsequent days. It was 
considered whether incorrect modelling of the reflected galactic emission could account for 
the discrepancy between predicted and observed Moon temperatures, however the spectral 
index required to fit the data is approximately f} — 4, which is unrealistically high. Given 
the clear presence of reflected terrestrial FM radio transmissions from the Moon, it is likely 
that this discrepancy is due to additional reflected transmissions from Earth, outside of the 
FM radio band. 



4. Analysis 

The strength of Earth's radio leakage signal reflected by the Moon is of interest for 
several reasons: 

(i) Reflected radio emission presents an additional challenge for experiments aiming to use 
the Moon as a smooth spectrum calibration source for detection of the EoR signal, 

(ii) The Moon will appear as a transient radio source that needs to be accounted for in 
current and future low frequency surveys and experiments, and 

(iii) The strength of the detected signal allows us to calculate how bright the Earth is as a 
radio source in this band and hence whether it might be detectable to technologically 
advanced extraterrestrial civilisations. 

To address point (iii) above, we can use these o bservations to calculate the equivalent 



isotropic power (EIP) of the Earth in the same way as lSullivan &: Knowlesl (119851 ) . However, 
because the MWA is an interferometer, our measurements only give us the difference between 
the background sky flux density and the Moon's flux density. The flux corresponding to the 
'zero spacing' is missing. We overcome this by converting the observed Moon flux density to 
brightness temperature by rearranging equation ([2]) and adding the known sky temperature 
given by equation (pQ). This gives us the observed Moon brightness temperature as shown in 
Figures O and [TUJ The brightness temperature of the Moon has been studied at a number o f 



wavelengths, however the most relevant to these observations is that of lHeiles &: Drake! (119631 ) 
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Fig. 7. — Integrated flux density of the 'Moon' pixels in the Moon-sky difference image on 
Day 1 when the Moon was at the centre of the beam and Day 2 when the Moon had drifted 
to approximately the half-power point of the beam. The flux density scale has been set using 
3C33 as a reference, but no primary beam correction has been applied to the Moon spectra. 



-17- 



3 

o 
c 

CD 


b 

'w 
c 

Q 

X 



5 - 



-10 



m r 



i — n r 



i 1 r 



i 1 r 




Observed 

Predicted 



J I L 



J I L 



J I L 



J I L 



80 



100 120 140 

Frequency (MHz) 



160 



Fig. 8. — Comparison between the observed Moon-sky flux density difference on 24 Septem- 
ber and the difference predicted for a model of the Moon that includes thermal emission 
and reflection of galactic background only. The excess flux density in the observations is 
attributed to the reflection of the Earth's radio leakage. 
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which shows that we should expect the Moon to have a constant brightness temperature of 
230 K at these frequencies. The excess brightness temperature due to reflected emission 
from the Earth is therefore calculated by subtracting the Moon's thermal emission and the 
component due to reflected galactic emission from this spectrum. This excess brightness 
temperature can then be converted back to flux density S m by equation (J2]) and to an EIP 
in Watts by: 

4 2 7rlO- 26 ^A^ m 

-L earth — 2 V / 

r mPm 

where d m is the average Earth-Moon distance in m, Au is the size of the frequency bins in 
Hz, S m is the observed flux density of th e Moon in J y, r m is the radius of the Moon in m, 
and p m is the Moon's reflectivity of 0.07 (jEvanslll969[ h 
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Fig. 9. — Disk average brightness temperature of the Moon calculated fr om MWA 32T mea 
surem ents, combined with a model of the sky temperature derived from iRogers fc Bowman 
J2008h and lde Oliveira-Costa eta3 fj2008h . 
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Fig. 10. — Disk average brightness temperature of the Moon calculated fr om MWA 32T mea- 
surem ents, combined with a model of the sky temperature derived from lRogers fc Bowman 
(120081 ) and Ide Oliveira-Costa et all ( 120081 ) . This is the same as Figure [9] but plotted on a 
different scale and with error bars based on the RMS noise of the MWA 32T difference 
image. The do t ted li ne is the thermal brightness temperature of the Moon of 230K from 



Heiles &: Drake! (11 9631 1 plus the expected reflection of galactic emission. 
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Evaluating equation (JSJ) for the Moon's spectrum using 5 MHz wide frequency bins gives 
the histogram in Figure [TTJ To understand what the signal shown in Figure [11] means from 
the perspective of the Search for Extraterrestrial Intelligence (SETI), we can calculate how 
bright the Earth might appear to an extraterrestrial observer. Summing the EIP over the 
20.5 MHz of the FM band gives a value of approximately 77 MW, so to an observer orbiting 
our closest neighbouring star, Proxima Centauri, the Earth would appear to be an extremely 
faint source of approximately 19 nJy within t he FM radio band. As a c omparison, the quiet 
Sun is a source of approximately 60,000 Jy ( lLantos fc Avignon! Il975l ) as viewed from the 
Earth in this frequency band during periods of low solar activity. The Sun's total EIP in 
the FM band is therefore approximately 3.5 GW, so viewed from Proxima Centauri it would 
appear as a 0.85 /xJy source and greatly outshine the Earth. 

Detection of such a small signal by a nearby, technologically advanced civilisation would 
require far superior instrumentation to our own, with even greater collecting area than is 
expected for future telescopes such as the Square Kilometre Array (SKA). For a three sigma 
detection at these frequencies with an SKA equivalent instrument with baselines long enough 
to resolve the Earth from the Sun, the integration time required is on the order of hundreds 
of years. This is due largely to the hig h background sky tempera ture which dominates the 
system temperature at low frequencies. lLoeb fc Zaldarriagal (120071 ) found that beamed radio 
leakage from a civilisation within 23 pc of Earth, using high power military type radars 
similar to our own, should be detectable by the next generation of low frequency radio 
interferometers. However we note that most of these types of systems on Earth operate at 
higher frequencies. For this reason, radio leakage from the Earth at higher frequencies is 
more likely to be detectable, as the background sky temperature is lower and our highly 
beamed, military radars operate in the GHz frequency range. SETI experiments p lanned 
for t he SKA are therefore more focussed at frequencies between 1 and 10 GHz ((Tarter 
2004 ). The long baselines of the SKA will also provide the required angular resolution and 
RFI rejection capabilities required for such SETI searches, the importance of which have 
been demonstrated by recent proof -of-concept SETI observa tions using Very Long Baseline 
Interferometry (VLBI) techniques (IRampadarath et alll2012l ). 



The observed Moon-sky flux density difference indicates that more observations are 
needed to properly characterise the Moon's radio properties at low frequencies if the Moon 
is to be used as a smooth spectrum calibrator for detection of the EoR signal. The Very 
High Frequency (V HF) part of th e radiofrequency spectrum from 30 to 300 MHz is used for 
many applications (lACMAl 120091 ) including FM broadcasting, aeronautical radionavigation 
and communications and amateur radio to name a few. These will all contribute to the 
reflected radio signal received on Earth from the Moon as evidenced by Figure [HJ However 
there are methods that could be employed to reduce the impact of this reflected RFI on 
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Fig. 11. — The equivalent isotropic power (EIP) of radio emissions from Earth as measured 
from MWA 32T observations of the Moon. The plot is a histogram with 5 MHz bins. 
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observations and some of these are discussed below. 

The RFI signals from the Moon are expected to be strongest from the sub-Earth point 
where the surface behaviour is closest to specular reflection. The higher spatial resolution of 
the completed MWA with 128 tiles and 3 km baselines will enable us to isolate the specular 
and diffuse reflection regions and take spectra from regions of the Moon less affected by 
reflected RFI from Earth. 

A higher spectral resolution investigation may also provide a means of av oiding reflected 



RFI. In many countries including Australia and the United States (|ITU||2001[ ) FM broadcast 
frequencies are assigned with a channel separation of 200 kHz and a maximum frequency 
deviation of ±75 kHz, leaving a 'guard band' of at least 25 kHz above and below the modu- 
lation. Observation in these guard bands should be free from RFI contamination or at least 
be affected at a reduced level. Similarly radio-quiet bands may be found in regions of the 
spectrum outside of the FM band if observing at higher spectral resolution. 

The reflected radio emission from the Moon will clearly also vary strongly as a function 
of time as different regions of the Earth come into the line-of-sight of the Moon. Since 
terrestrial broadcasting is generally directed toward the horizon rather than straight up, the 
worst RFI contamination is likely to occur when the Moon's position falls near the plane of 
the horizon at the most highly populated regions of the Earth. Hence RFI contamination 
may be reduced by planning observations at times when the geometry between the Moon, 
the telescope and population centres on Earth minimises the radio signal beamed toward the 
Moon. 

An additional concern is that the occasional presence of the Moon in observations of 
deep fields that are visited repeatedly may present complications due to the rapid movement 
of the Moon with its strong reflected signal, which itself is rich in time-variable spectral 
structure as the Earth rotates and different population centres illuminate the Moon with 
radio transmissions. Furthermore, the Moon subtends sufficient solid angle that ensembles 
of background sources will be steadily occulted and re-exposed in areas of long integration. 
This is a concern for some instruments, since the principles of aperture synthesis and co- 
herent integration with interferometers rely on the constancy of the radio sky throughout 



the in tegration. However, instruments such as the MWA and LEDA (jGreenhill fc Bernardi 



20121 ) will implement snapshot imaging with subsequent combination in the image domain, 



an approach that might be made robust to this kind of variability. 
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5. Conclusions 

The MWA 32 tile prototype system was used to observe the Moon over two subsequent 
nights in order to test the viability of using the Moon as a smooth spectrum calibration 
source for EoR signature detection. A direct differencing experiment between the Moon and 
sky flux density from the same patch of sky over the two nights found that the Moon does 
exhibit the properties expected for a cool thermally emitting body, but that the spectrum is 
corrupted by radio leakage from Earth, most notably within the FM radio band of 87.5 to 
108 MHz. From this reflected signal we were able to measure the EIP of the Earth in the 
FM radio band which was found to be approximately 77 MW. 

The implications of these observations are that RFI from Earth will be a significant 
challenge to overcome if the Moon is to be used as a smooth spectrum calibration source 
for EoR detection. Reflected terrestrial radio emissions will have to be carefully taken into 
account in any future experiments and more Moon observations are required to characterise 
this RFI. Additionally, it was found that the total signal being transmitted from the Earth 
in the FM radio band would be difficult to detect, even by a nearby technologically advanced 
civilisation with superior radio astronomy technology to our own. 
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